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1.0  INTRODUCTION 

This  report  presents  theoretical  techniques  which  may  be  used  in  the 
analysis  and  design  of  microstrip  antennas.  The  involvement  of  the  Physical 
Science  Laboratory  in  the  theoretical  analysis  of  these  electrically  thin 
radiators  has  been  principally  motivated  by  a practical  requirement  for 
better  design  techniques  which  could  be  applied  to  the  development  of 
microstrip  antennas  for  systems  applications.  Over  the  past  seven  years, 

PSL  has  acquired  considerable  practical  experience  in  the  design  of  both 
stripline  and  microstrip  antennas;  microstrip  antennas  have  been  produced 
in  wrap-around  cylindrical  arrays  for  sounding  rocket  application 
[Weinschel,  1975]  and  in  planar-array  versions  for  UHF  communication  from 
buoys  and  ships  to  geosynchronous  satellites  [Weinschel  and  Carver,  1976]. 
Until  recently,  however,  these  designs  have  come  from  a combination  of 
experimental  cut-and-try  techniques  supplemented  by  standard  transmission- 
line theory  applied  to  simple  rectangular  patches. 

The  microstrip  antenna  concept  dates  back  some  twenty-five  years  to 
work  by  Deschamps  [1953].  However,  nearly  eighteen  years  passed  before  the 
work  was  renewed  and  reported  by  J.  Q.  Howell  at  NASA  Langley  Research  Center 
[1972,  1975],  R.  E.  Munson  at  Ball  Brothers  Research  Corporation  [1972, 

1975]  and  H.  D.  Weinschel  at  the  Physical  Science  Laboratory  [1973,  1975]. 

Mathematical  modeling  of  the  basic  microstrip  radiator  was  initially 
carried  out  by  applying  transmission-line  analogies  to  simple  rectangular 
patches  fed  at  the  center  of  a radiating  wall  [Munson,  1974;  Derneryd,  1976]. 
The  radiation  pattern  of  a circular  patch  was  analyzed  and  measurements 
reported  by  Carver  [1976]. 

These  earlier  mathematical  models,  while  giving  reasonably  good  results 
for  center-fed  rectangular  patches,  all  shared  some  serious  deficiencies: 

(1)  they  consistently  predicted  resonant  frequencies  which  were  a few  percent 
off  (usually  on  the  high  side)  from  measured  resonant  frequencies,  (2)  they 
could  not  be  adapted  to  multi-mode  structures  with  arbitrary  feed  points 
(such  as  corner-fed  rectangular  patches),  and  (3)  they  could  not  be  adapted 
to  patches  of  arbitrary  shape. 

The  present  analytical  work  began  in  January,  1977  and  received  ARO 
support  beginning  in  April,  1978.  It  has  been  directed  along  two  comple- 
mentary paths:  (a)  for  rectangular  and  circular  patches,  a modal  expansion 
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method  incorporating  impedance  boundary  conditions  at  the  radiating  walls, 
and  (b)  for  patches  of  arbitrary  shape,  a numerical  approach  based  on  the 
method  of  finite  elements.  The  initial  proposal  for  the  numerical  analysis 
was  to  use  the  Unimoment-Monte  Carlo  technique  [Coffey,  1976;  Coffey  and 
Carver,  1977];  however,  for  reasons  discussed  in  Sec.  3.0,  this  numerical 
technique  has  been  abandoned  in  favor  of  the  method  of  finite  elements. 

During  the  course  of  this  effort,  a comprehensive  report  was  received  from 
workers  at  the  University  of  Illinois  [Lo,  et.  al. , 1977]  which  used  a 
similar  expansion  method  applied  to  impedance  and  radiation  pattern  analysis 
of  rectangular,  circular  and  triangular  patches.  One  of  the  differences 
between  the  Illinois  report  and  the  present  approach  is  that  the  latter 
incorporates  impedance  boundary  conditions  at  the  radiating  walls,  whereas 
the  former  does  not.  The  wall  conductance  used  in  the  present  report  is 
based  on  a parallel-plate  TEM  waveguide  radiating  into  a half-space 
[Marcuvitz,  1949;  Harrington,  1961].  The  wall  susceptance  is  based  on  a 
formulation  by  Hammerstad  [1975].  The  result  of  this  model  is  a prediction 
of  impedance  vs.  frequency  behavior  which  shows  very  good  agreement  with 
measured  results  for  rectangular  patches;  there  is  also  good  agreement  with 
measured  results  for  circular  patches. 

This  report  begins  with  a discussion  of  a general  mathematical  frame- 
work for  the  analysis  of  microstrip  antennas.  A discussion  is  then  presented 
of  methods  based  on  the  finite  elements  approach  along  with  key  aspects  of 
a computer  code  developed  along  these  lines.  Next,  a classical  modal 
expansion  method  is  presented  and  applied  to  both  rectangular  and  circular 
patches.  Finally,  the  finite-elements  approach  is  used  to  solve  for  the 
magnetic  currents  on  the  radiating  walls  of  a five-sided  patch  which  can 
produce  circular  polarization. 
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2.0  THEORETICAL  MODELS 

2 . 1 Introduction 

A microstrip  antenna  is  shown  in  Figure  2.1  as  a patch  of  arbitrary  shape 
separated  from  a large  ground  plane  by  a homogeneous  dielectric  of  thickness 
t,  where  t is  assumed  to  be  much  less  than  a wavelength.  Practical  microstrip 
antennas  have  been  fabricated  with  rectangular,  circular,  semi-circular,  tri- 
angular and  pentagonal  patch  shapes.  It  is  therefore  of  interest  to  develop  a 
theoretical  foundation  which  can  handle  patches  of  arbitrary  size  and  shape. 

The  objectives  of  such  an  analysis  are  to  predict  both  the  far-field  pattern 
characteristics  (including  gain  and  polarization)  as  well  as  such  near-field 
characteristics  as  input  impedance,  impedance  bandwidth,  antenna  efficiency, 
mutual  coupling,  etc.  This  is  an  ambitious  program  and  will  require  the 
efforts  of  many  investigators  before  it  is  completed.  It  is  expected  that 
approximation  methods  will  be  needed  in  order  to  obtain  practical  answers . 
However,  to  date  there  have  been  no  publications  in  the  open  literature  which 
attempted  to  set  up  a formalism  for  the  microstrip  antenna  based  on  a rigorous 
approach.  It  seems  worthwhile  to  attempt  such  an  attack  even  though  satisfac- 
tory solutions  may  be  difficult. 

The  analysis  of  the  microstrip  antenna  may  be  facilitated  by  dividing  the 
radiation  into  an  interior  region  and  an  exterior  region.  The  interior  region 
is  a cavity  formed  by  the  patch,  the  portion  of  the  ground  plane  under  the 
patch,  and  the  radiating  walls  formed  by  the  projection  of  the  patch  perimeter 
onto  the  ground  plane.  The  exterior  region  is  the  rest  of  space  and  includes 
the  remainder  of  the  ground  plane,  the  remainder  of  the  dielectric  and  the  top 
of  the  patch  conducting  surface.  Practical  ground  planes  are  usually  several 
wavelengths  in  dimension  so  that  the  assumption  of  an  infinite  ground  plane  is 
justified  for  purposes  of  computing  the  far-field  pattern  shape  over  most  of 
the  upper  hemisphere.  However,  for  angles  near  the  plane  of  the  ground  plane, 
edge  diffraction  effects  must  be  considered. 

The  connection  between  the  interior  and  exterior  problems  may  be  made  by 
enforcing  the  boundary  condition  that  there  be  continuity  of  the  tangential 
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components  of  E and  H at  the  radiating  walls  of  the  cavity.  This  is  an  extremely 
difficult  problem  and  nas  not  to  date  been  solved. 

2 . 2 Notation 

A notation  is  needed  which  explicitly  identifies  both  interior  and  exterior 
fields.  If  the  distinction  between  interior  and  exterior  fields  is  not  obvious 
from  context,  the  symbol  i used  as  a superscript  will  identify  an  interior 

— — i — 0 -g 

field,  i.e. , E or  H . Similarly,  (E  , H ) are  exterior  fields. 


S I 


Subscripts  are  then  used  to  denote  either  mode  indices  (vide  Section  3.0)  or 

the  value  of  a field  at  a particular  location.  In  particular,  the  fields 

evaluated  in  the  thin  ribbon  forming  the  radiating  walls  are  denoted  as  (E  , 

H ).  Fields  evaluated  on  the  exterior  ground  plane  are  denoted  as  (E  , H ) 
rw  6 gp  gp 

and  exterior  fields  evaluated  on  the  top  of  the  radiating  patch  are  denoted  as 

(E  , H ). 

P P 

As  in  most  other  published  literature,  primed  coordinates  are  used  to  denote 
source  point  location  and  unprimed  coordinates  are  used  to  denote  field  point 
location. 

2 . 3 Formalism  for  the  Interior  Problem 

Consider  the  interior  region  of  a cavity  with  radiating  side  walls,  as  shown 
in  Figure  2.2.  Since  t « A.,  it  may  be  assumed  that  the  interior  electric 
field  is  entirely  z-directed  and  is  independent  of  z. 

The  interior  fields  must  satisfy  the  wave  equation 


V x V x E - k2E1  = - jui|j  f - V x Hs 


V x V x H1  - k2H1  = - ju)£  MS  + V x Is 


(2.1a) 


(2.1b) 


where 


k2  = ui2  u e 
ro 


and 


'C' 


e = £'  - js"  = £Q(£j.  - j e'p  . (2.2) 

— S — 2 

J is  the  electric  current  density  (Am  ) due  to  the  coaxial  feedline  source 
and  may  be  written  as 

js  = z — ■ 6(x  - xq)  6(y  - yQ)  (o  < 2'  < t)  (2.3) 

a 


where  a is  the  cross-sectional  area  of  the  coaxial  current  element  which  is 

-s  -2 

assumed  to  be  electrically  thin.  M is  the  magnetic  current  density  (V  m ) 

due  to  the  microstrip  feedline  source  and  may  be  written  as 


MS  = < 


v(*' >0)  a 


o(y) 


tw 


f W , Wx 

(x < x'  < X + -) 

2 2 

(o  < z'  < t) 
otherwise 


(2.4) 


where  w is  the  width  of  the  microstrip  feedline  and  V(x' ;o)  is  the  voltage 
between  the  patch  and  the  ground  plane  evaluated  at  the  feedline  connection 
point. 


For  an  antenna  fed  only  by  a microstrip  feedline,  J = 0;  similarly,  for 

—5 

an  antenna  fed  by  a coaxial  feeline,  M =0. 


In  addition  to  satisfying  the  wave  equations  (2.1),  the  fields  must  also 
satisfy  the  following  boundary  conditions*: 


7TL 

n x E =0 


(on  top  and  bottom  conductors) 


(2.5) 


rre  _ Tre 

n x E = n x E 

rw  rw 


771  re 

n x H = n x E 
rw  rw 


(on  radiating  walls) 


(2.6) 


*This  first  boundary  condition  assumes  that  the  patch  and  the  ground  plane  are 
perfect  conductors.  Since  almost  all  practical  microstrip  antennas  are  fabricated 
from  copper,  this  is  an  excellent  assumption. 
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If  the  interior  region  were  a non-radiating  cavity,  then  H = o so  that  the 
tangential  components  of  H1  at  the  radiating  walls  would  be  zero.  In  the 
language  of  classical  boundary-value  problems,  this  is  known  as  a Neumann  con- 
dition. 

However,  the  practical  problem  must  contend  with  radiating  walls  so  that  the 
interior  fields  must  satisfy  an  impedance  condition.  A solution  for  the  fields 
then  requires  a knowledge  of  the  impedance  presented  by  the  radiating  walls  as 
a function  of  location  on  the  wall. 

The  classical  method  for  handling  such  cavities  is  to  find  a complete  set  of 
eigenvectors  which  can  be  used  to  expand  the  general  field.  This  is  not  diffi- 
cult for  such  simple  geometries  as  rectangular  and  circular  patches  and  has 
been  discussed  at  some  length  by  Lo,  et.  al.  [1977].  If  a complete  basis  can- 
not be  found,  approximation  techniques  such  as  the  Ritz-Galerkin  method,  point 
matching,  method  of  moments,  etc.,  can  be  used  to  find  the  interior  fields. 

Section  3 of  this  report  describes  an  approach  to  the  interior  problem  using 
the  method  of  finite  elements.  This  method,  which  shows  great  promise,  requires 
a knowledge  of  the  impedance  of  the  radiating  walls  or  a reasonable  approxima- 
tion to  it. 

The  balance  of  complex  power  for  the  interior  region  may  be  described  by 


P.  = PD  + PJ  + 2jw  (W  - W ) 
in  R d J m e 


(2.7) 


where 


P^n  = complex  power  supplied  by  source 

P^  = complex  power  radiated  through  edge  walls 

Pj  = power  dissipated  by  conducting  walls  and  by  the  dielectric 

w = angular  radian  frequency 

Wm  = time-average  magnetic  energy  stored  in  cavity 
W = time-average  electric  energy  stored  in  cavity 


2.3.1  Complex  Source  Power 


For  microstrip  antennas  fed  by  coaxial  lines,  the  input  power  may  be  computed 
by 


P" 


- = ' Iff  r 


J dv' 


(2.8) 


Using  (2.3)  for  J and  assuming  E1  = z E1,  (2.8)  becomes 


v 

^c)  = - / EL(x  , y ) IW 

in  /o’  Jo 


(z ' )dz ' 


(2.9) 


If  I(z')  is  a constant,  then  (2.9)  becomes 


= v.  ii 

in  in  in 


(2.10) 


where 


V.  = 
in 


J £l  (x0’  yo)dz’  = ' c £l(xo’  yo} 


(2.11) 


Because  of  the  assumption  that  J was  filamentary  (c.f.  eqn.  2.3)  the  reactive 
input  power  computed  as  the  imaginary  part  of  (2.9)  will  be  in  error.  This 
will  be  discussed  later. 


For  microstrip  antennas  fed  by  microstrip  feedlines,  the  input  power  may  be 
computed  by 


.(ms) 


in 


= .fff 


(2.12) 


,(ms)  _ 


in 


■I 


and  assuming  H1 

= x + y »;> 

+ w/2  t 

(H1)"'  vte.'.i.°)  dx-  dz< 

(2.13) 

- w/2  *4 

x tw 

[ 

. - 

>*«•* 

- - - - - 


~yr 
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If  the  strip  width  is  narrow,  V(x' ;o)  is  essentially  constant  and  equal  to 
V^in  so  that  (2.10)  becomes 


p(ms)  _ j*  v 

in  in  in 


(2.14) 


where 


(2.15) 


2.3.2  Dissipated  Power 

The  power  dissipated  as  heat  can  be  subdivided  into  that  associated  with  con- 
ductor losses  (P.  ) and  that  associated  with  dielectric  losses  (P,.) 

ac  dd 


P . = P , + P,  , 

d dc  dd 


(2.16) 


where  [Kraus  and  Carver,  1973]  : 


P . = Re  Z 

dc  c 


ds 


(2.17) 


cond . 
walls 


and 


Pdd  = W£ 


" ///'Ei| 


dv 


(2.18) 


cavity 


2.3.3  Radiated  Power 


The  complex  radiated  power  is  obtained  by  integrating  the  complex  Poynting 
vector  over  the  radiating  wall  surface,  i.e. 


jj  (E1  x H1*) 


ds 


(2.19) 


rad. 

wall 
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Since  E1  is  entirely  z-directed,  the  only  component  of  H1  which  is  of  interest 
is  that  tangent  to  the  radiating  wall  and  parallel  to  the  local  edge  of  the 
patch;  this  component  is  designated  as  H^.  Thus, 

h = - f vt  K dz  <2-2o> 

wall 

perimeter 

where  is  the  voltage  across  the  radiating  wall  and  is  the  parallel  com- 
ponent of  the  magnetic  field.  is  also  related  to  the  surface  electric 

current  density  at  the  edge,  as  shown  in  Figure  2.3,  by  the  boundary  condition 

J = n x H„  (2.21) 

where  n is  a unit  nor.ial  directed  away  from  the  conducting  surface.  If  is 
entirely  tangent  to  tne  radiating  wall,  then  the  current  on  the  bottom  side  of 
the  patch  will  flow  around  the  edge  and  onto  the  top  side  of  the  patch,  as 
shown.  The  electric  field  just  exterior  to  the  radiating  wall  will  fringe  and 
refract  at  the  dielectric-air  interface  in  accordance  with  Snell's  Law. 


2.3.4  Stored  Energy 


The  total  stored  energy  in  the  cavity  is  given  by 
W = Wm  + We  = “ jff  (£Ie1!  2 + M0 1 H1!  2,dv 


cavity 


(2.22) 


Since  the  interior  fields  are  independent  of  z,  this  can  be  slightly  simpli- 
fied to 


W 

e 


// 

cavity 

area 


iE^I2  ds 


W 

m 


l Mo  C JJ  IhV  ds 

cavity 

area 


(2.23) 


(2.24) 
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2.3.5  Input  Impedance 


The  input  impedance  to  a coaxially-fed  microstrip 

2 

computed  by  noting  the  = j I j Z and  using  (2-9) 


antenna  can  be 
to  obtain 


Z 


in 


A. 

s 


Ei(x 


,y0)i*(z') 


dz ' 


(2-25) 


where  a z-directed  filamentary  probe  has  been  assumed.  Similarly, 

the  input  admittance  to  a microstrip  transmission  line-fed  microstrip 

2 

antenna  can  be  computed  by  noting  that  P.  = |V.  I Y.  and  using  (2-12) 

in  \ m l in 

to  obtain 


Y. 

in 


[Hi(x1,y1)]2  V(D  d SL 


(2-26) 


where  the  feedline  input  has  been  modeled  as  a thin  strip  of  magnetic 

“■  A A /\ 

surface  current  density  Mg  = (z  x n)(V/t)  where  n is  a unit  vector 
normal  to  the  thin  radiating  surface  and  pointing  into  the  cavity.  The 
strip  is  of  width  w.  If  the  strip  width  is  narrow,  then  V(£)  is  essen- 
tially constant  and  equal  to  V . If  t is  small,  then  the  input  current 
on  the  lower  side  of  the  upper  microstrip  transmission-line  conductor 
is  Iin  = [H1(x  ,y  )]w,  so  that 


in 


(2-27) 


Likewise,  from  (2-25),  if  a small  value  of  t is  assumed  then  both  E and 
I(z')  are  constant  so  that 


'in 


v,  /I, 
in  in 


(2-28) 


The  usual  procedure  for  a probe-fed  cavity  is  to  find  V « -tE^x^.y^ 
by  expanding  Ez(x,y)  in  terms  of  a complete  set  of  orthogonal  modal 
vectors  (eigenfunctions)  emn(x»y)  an^  using  the  orthogonality  property  to 
find  the  coefficients  of  the  expansion  in  terms  of  the  input  current. 

A similar  procedure  can  be  used  for  the  microstrip-line  fed  cavity  by 
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finding  I = [H^Cx^.y^lw  and  expanding  H*. 

In  this  classical  eigenvalue  problem,  a denumerably  infinite  set 

of  discrete  eigenvectors  and  eigenvalues  are  found  by  imposition  of 

either  Dirichlet  (E  = 0 on  boundaries)  or  Neumann  (^E  fan  = 0 on 
z z 

boundaries)  boundary  conditions.  The  mode  vectors  are  customarily 
orthonormalized  as  follows: 


cavity 

volume 


5 

mn,  pq 


where  & 


mn,pq 


1 if  m = p and  n = q;  <T  = 0 otherwise.  Thus, 


(2-29) 


»* 

I 

i 

I 

I 

I 


dv  ■ 1 (2-30> 

t h 

is  the  electric  stored  energy  in  the  mn  mode.  In  a closed  cavity,  the 

time-average  electric  and  magnetic  energies  are  equal.  Therefore, 

normalization  of  the  e^  vectors  automatically  normalizes  the  corresponding 

h vectors, 
mn 

When  the  cavity  begins  to  radiate  real  and  reactive  power  out  the 
radiating  walls,  several  differences  are  noted: 


1.  The  time-average  electric  and  magnetic  energies  are  no  longer 
equal . 

2.  The  eigenvalues  are  no  longer  necessarily  real  or  discrete. 

3.  Energy  is  cross-coupled  between  discrete  modes. 

4.  Only  a finite  number  of  discrete  eigenfunctions  exist. 

5.  In  addition  to  these  discrete  eigenfunctions,  there  will  be 
added  a continuously  varying  eigenfunction  with  a continuous 
range  of  eigenvalues. 

The  electric  field  Ez(x,y)  interior  to  a probe-fed  cavity  satisfies 
the  inhomogeneous  differential  equation: 


A A 


k2  E 


iy*b  xin  * (x  " xo)  ^ (y  ’ yo}  (2'31) 
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* 
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F 


■ 


Tfc’ 


along  with  the  impedance-type  boundary  conditions: 


2>EZ 

+ oC  E » 0 (on  radiating  walls)  (2-32) 

2>n 


The  partial  derivative  S/3  n is  with  respect  to  an  outward-directed 

coordinate  and  oC  is  dependent  on  the  impedance  of  the  walls.  The 

complete  set  of  functions  required  to  represent  any  arbitrary  E^ 

consists  now  of  both  the  discrete  modes  e (x,y)  plus  a continuous 

mn 

spectrum  £(x,y;k  ,k  ).  Thus  the  solution  is  given  by 
x y 


Ez(x,y) 


-JOJfo  Tin  ^ G(x,y;x'  ,y')S  (x'-xq)£  (y'-yQ)  dx'  dy' 


cavity 


(2-33) 


or 


Ez(x,y)  - -ju iy.o  Iin  G(x,y;xo,yo) 


(2-34) 


where  G is  the  Green's  function  which  satisfies  the  relation 


i.U 


M N 


G(x,y|x  ,y  ;k  ,k  ) dk  dk  = 5*  e (x,y)e  (x  ,y  ) 

o o x y x y Cj„  mn  J mn  o J o 


C1’C2 


m=0  n=0 


IS  t<x,y;k-' 


;k,k)e(x,y;k,k)dk  dk 
x y o Jo  x y x y 

(2-35) 


where  the  contour  over  k is  determined  by  that  over  k according  to 

y 22  1/2  x 

(k  - k__)  . The  contour  (or  C2)  is 


the  separation  condition  k 

y x 

taken  around  all  singularities  of  the  Green's  function.  Equations  (2-31) 
through  (2-35)  are  developed  in  the  context  of  a rectangular  microstrip 
antenna;  a similar  development  for  circular  patches  can  be  carried  out 
using  cylindrical  coordinates.  For  both  of  these  cases,  the  two- 
dimensional  Green's  function  is  separable  into  a product  of  one-dimensional 
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Green's  functions.  However,  for  arbitrary  microstrip  patch  shapes, 
this  is  not  in  general  true. 

Returning  to  (2-35),  it  is  noted  that  for  closed  cavities,  the 
complete  set  of  mode  vectors  is  represented  by  only  the  discrete  sum. 
For  open  radiating  cavities  (waveguides)  which  satisfy  a radiation 
condition  at  the  walls,  the  complete  set  of  modes  is  the  continuous 
eigenfunction  seen  in  the  contour  integral  on  the  right-hand  side  of 
(2-35).  For  the  most  general  case,  however,  both  discrete  modes  and 
continuous  modes  will  be  present,  along  with  discrete  and  continuous 
eigenvalues. 

Practical  microstrip  antennas  usually  have  impedance  bandwidths 
which  are  small  and  of  the  order  of  1%.  This  indicates  a high-Q 
condition  and  that  the  interior  fields  are  dominated  by  the  discrete 
modes.  However,  the  continuous  spectrum  (“(x.yjk^.k^)  cannot  be 
ignored  since  it  is  closely  related  to  the  radition  admittance  of  the 
walls. 

These  impedance  considerations  will  be  applied  to  specific 
microstrip  patch  shapes  in  Section  4 of  this  report. 

2.3.6  Impedance  Bandwidth 


The  impedance  bandwidth  is  given  by 
f 

Af  = — 


(2-36) 


where  fr  is  the  resonant  frequency  and  the  total  Q is  given  by 
W 


(2-37) 


Re  P 


in 


W + W being  the  total  stored  energy  both  in  the  cavity  interior 
e m 


with  W 

and  the  near-field  region  of  the  antenna.  The  real  input  power  is  either 
dissipated  by  the  conductors  or  dielectric  or  is  radiated  as  real  power 
through  the  radiating  walls.  Thus 


(2-38) 


where 


Qr  = “ RTF 


(2-39) 


Q =oj  

Re  Z // | H |2  ds 
c J J tan  1 


(2-40) 


d e^/l/lE1!2  dv 


(2-41) 


The  total  stored  energy  (both  in  the  cavity  and  in  the  radiated  field) 
can  by  computed  by 


W = y ///(elE1!2  + u | H1 1 2)  dv  + ^ Im  Pr 


(2-42) 


where  P^.  is  given  by  (2-20) 


2.3.7  Antenna  Efficiency 


The  antenna  efficiency  is  the  ratio  of  real  radiated  power  to  real 
input  power,  i.e. 


Thus, 


(2-43) 


n = Q/Q, 


(2-44) 


Also,  using  (2-20),  (2-43)  can  be  written  as 


R (R  2 + X 2)  1 

n = — is iS_  x x | V.  | 2 dH 

R,  (R2  + X2)  t |V.  |2 
in  w w 'in' 


(2-45) 


where  V.  = -E  t evaluated  at  the  radiating  wall,  and  where  Z = R + jX 

l Z W W J W 

is  the  wall  impedance,  assumed  independent  of  perimeter  position. 
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2.4  Formalism  for  the  Exterior  Problem 

The  two  objectives  of  the  exterior  problem  are  to  find  the  far-field 
radiation  pattern  (including  gain  and  polarization)  and  to  describe  the 
near-field  coupling  behavior  of  the  microstrip  patch.  The  exterior  fields 
may  be  found  from  the  equivalence  principle  if  the  tangential  fields  in 
the  radiating  walls  of  the  cavity  and/or  the  tangential  magnetic  fields 
on  the  exterior  conductors  are  known;  the  equivalent  magnetic  surface 
currents  in  the  gap  are  given  by 

Mg  = E x n (2-46) 

where  is  evaluated  in  the  gap  (radiating  wall)  and  n is  a unit  vector 

which  is  normal  outward  from  the  gap.  The  equivalent  electric  surface 
currents  on  the  top  of  the  patch  and  on  the  ground  plane  are  given  by 

J = z x h (2-47) 

s 

where  H is  evaluated  on  either  the  patch  or  the  ground  plane. 

The  exterior  fields  are  then  found  by  the  formalism 

E6  = ~ — (VxVxA  - J)  - Vx?  (2-48) 

jwe 


When  the  exact  current  densities  are  known  and  if  the  conductivities  of 
the  patch  and  ground  plane  are  infinite,  the  fields  satisfy  the  boundary 
conditions: 


z 


0 


on  patch  and  on  ground  plane 


(2-52) 


x E6 
z * 

n x H6 

e*  x ; 

lim  |rEe|  is  finite 
r-*“ 

radiation  condition 

lim  |rHe|  is  finite 

r-x= 


= J 


~ — i 

= n x H 


on  radiating  wall 


= M 


s J 


(2-53) 


(2-54) 


The  formalism  of  (2-48)  - (2-51)  requires  a knowledge  of  both  Jg 

and  M , as  shown  in  Fig.  2.4a.  The  fields  can  also  be  found  from  J 
s _ s 

alone  (Fig.  2.4c)  or  Mg  alone  (Fig.  2.4b). 

For  the  Mg  formulation,  a perfect  conductor  is  placed  over  the  surface 
bounding  the  interior  region.  Since  a tangential  Jg  over  a perfect 
conductor  produces  a zero  net  field  in  the  exterior  region,  the  field  is 
calculated  from  Mg  alone  in  the  presence  of  the  conductor.  The  exterior 
field  is  then  found  as 


E*(r) 


rad. 

wall 


Ez(r’) (z  x A) 
4r  |r  - r’ | 


ds  1 


(2-50) 


subject  to  the  boundary  condition 


n x E 


on  perfect  conductors 


(2-51) 


This  is  an  awkward  equation  to  solve  because  of  blockage  introduced 
by  the  conducting  step  jump. 

In  Fig.  2.5,  additional  equivalent  sources  are  shown  which  produce 
the  same  field  in  the  half-space  above  the  plane  of  the  patch  and  a null 
field  below  this  plane.  Fig.  2.5c  shows  an  equivalence  which  takes  advan- 
tage of  image  theory;  in  this  case,  the  perfect  conductor  is  removed  and 
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PERFECT  MAGNETIC  CONDUCTOR 


(c) 


Fig.  2.4.  Three  source  formulations  which  produce  the  same  exterior 
field.  (a)  J and  M ; (b)  M alone  plus  perfect  electric 
conductor;  (cf  J along  plusSperfect  magnetic  conductor. 
The  exterior  region  includes  the  half-space  above  the 
patch  as  well  as  the  region  interior  to  the  dielectric. 
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NULL  FIELD 
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Fig.  2.5.  Three  additional  equivalent  source  formulations  which 
produce  the  same  field  above  the  plane  of  the  patch, 
(a)  J and  M ; (b)  M alone  plus  perfect  conductor; 
(c)  2fl  alonl  (magnetic  image). 


the  magnetic  surface  current  density  is  doubled.  Thus, 


E*(r) 


S-S 

o 


(E  *z)e-jk|r-r 
§2 

2tt  |7  - r'  | 


ds' 


(2-52) 


where  the  aperture  plane  is  S and  the  portion  of  S which  represents 

the  patch  is  S . 

o 

Finding  an  exact  relationship  between  the  tangential  electric  fields 
in  the  aperture  plane  and  those  in  the  radiating  wall  of  the  cavity  is  not 
a simple  task,  as  is  illustrated  in  Fig.  2.6  where  the  electric  field  lines 
have  been  determined  by  graphical  field  mapping  procedures  for  the  d.c.  case. 
In  the  immediate  vicinity  of  the  edge,  these  d.c.  field  lines  will  be 
similar  to  those  of  a radiating  cavity.  Due  to  the  dielectric/air 
discontinuity,  the  field  lines  will  refract  according  to  Snell's  Law 
and  the  electric  field  in  air  will  be  larger  than  that  of  the  same 
line  in  the  dielectric.  Furthermore,  the  tangential  component  of  E 

ap 

is  appreciable  only  very  near  the  edge  of  the  patch. 

An  approximate  relation  between  the  field  z E in  the  radiating  wall 

rw 

of  the  cavity  and  the  tangential  field  E in  the  aperture  plane  may  be 
found  by  first  assuming  that  the  region  of  the  aperture  plane  over  which 
effective  magnetic  currents  exist  is  a thin  strip  of  width  w extending 
from  the  edge  of  the  patch  outward.  Over  this  region,  z-components  of 
E tend  to  cancel,  leaving  only  a tangential  component.  Let  be  the 
length  of  a field  line  from  the  ground  plane  to  the  dielectric  surface 
and  let  be  the  remaining  length  of  the  line  from  the  dielectric  to  the 
patch.  Then  assuming  the  patch  to  be  an  equipotential  surface. 


V = -(E,  I,  + E l)  (2-53) 

P Q Q 3.  3 

where  E,  is  the  electric  field  in  the  dielectric  and  E is  the  electric 
d a 

field  in  air,  both  evaluated  on  the  same  line.  Furthermore,  the  electric 
flux  density  is  constant  over  any  one  line  so  that 

D = eeE«eE  (2-54) 

o r d o a 
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GROUND  PLANE 


RADIATING  WALL 
OF  CAVITY 


Fig.  2.6.  Illustrating  the  d.c.  fringing  field 

distribution  for  the  edge  of  a microstrip 
patch  with  e = 2.6,  t = 0.0625"  and  a 
patch  conductor  thickness  of  0.002". 
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where  er  is  the  dielectric  constant  of  the  patch.  Combining  (2-54)  and 
(2-53), 


I 

I 

I 

I 

I 

I 
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v ( — + t r1 

p c a 

K r 


(2-55) 


At  the  edge  of  the  patch,  t is  zero  and  E is  entirely  tangential  so  that 

a a 


V 

E (edge)  = cr  f 
d 


e E 
r rw 


(2-56) 


This  is  the  maximum  tangential  field  in  the  aperture  plane.  The  effective 

tangential  field  is  less  than  this.  For  an  aperture  plane  of  effective 

strip  width  w = t,  the  tangential  field  at  the  outer  edge  of  the  strip 
• ★ 

is  essentially  zero  , so  that  the  effective  aperture  tangential  field  is 
approximately  half  that  given  by  (2-56) : 

_ er 

(over  aperture  strip  of  width  t)  (2-57) 


ap 


rw 


where  n^  is  a unit  vector  lying  in  the  aperture  plane  and  extending  out- 
ward from  the  patch  edge. 

The  far  fields  may  now  be  computed  by  letting  | r — r 1 | = r in  the 
denominator  of  (2-52)  and  by  letting  |r  - r'|  = r - r'  cos  <J>'  in  the 
exponential  of  (2-52)  ; these  coordinates  are  shown  in  Fig.  2.7.  Then, 
the  far  electric  fields  are  given  by 


E(9,<j>) 


9 E0(9,4>)  + 


E , (9,  ) 


where 


E9  * F» 


jk  F. 


(2-58) 


(2-59) 


— e 
c “*  — 


-jkr 


2t rr 


M_(r\*')  e^r'cos*'  i£, 


perimeter 


and  where 


(2-60) 


M = 


(z  x n ) 
P 


rw 


(2-61) 


It  is  emphasized  that  the  magnetic  surface  current  density  in 
(2-61)  is  approximate;  more  work  is  required  to  obtain  better 
relationships  between  the  magnetic  surface  current  density  in  the 
aperture  plane  and  the  z-directed  field  in  the  radiating  wall  of  the 
cavity. 

The  corresponding  far-zone  magnetic  fields  are  given  by 


H , - — 


H„ 


(2-62) 


where  Zq  = 376  0 is  the  free-space  impedance. 

Finally,  it  is  observed  that  the  real  radiated  power  obtained  by 
integrating  the  real  part  of  the  power  density  around  the  radiating  wall 
of  the  cavity  (eqn.  2-20)  will  be  equal  to  the  real  radiated  power 
obtained  by  integrating  the  Poynting  vector  over  a far-field  hemisphere 
minus  the  real  power  appearing  in  the  surface  wave  traveling  in  the 
dielectric . 
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3.0  NUMERICAL  TECHNIQUES  FOR  THE  INTERIOR  PROBLEM  • 

It  was  proposed  originally  that  the  most  efficient  numerical  method 
for  solving  the  interior  problem  was  the  Unimoment-Monte  Carlo  Method 
[Coffey,  1976;  Coffey  & Carver,  1977]  implemented  on  a medium-size 
hybrid  computer.  This  technique  was  tested  during  the  period  10  April, 

1978  - 30  June,  1978  and  was  found  to  perform  quite  well  on  simple 
boundary  value  problems  [Coffey  & Carver,  1978].  However,  it  has  been 
decided  to  abandon  this  unique  method  in  favor  of  more  traditional 
techniques  because  (1)  the  hybrid  computer  at  NMSU  was  being  repaired 
during  the  period  30  June,  1978  - 15  September,  1978  and  was 

unavailable  for  research,  and  (2)  one  of  the  principal  investigators  i 

(ELC)  terminated  employment  with  NMSU  and  a more  traditional  technique 

seemed  to  be  required  which  would  not  depend  on  his  expertise  in  the 

hybrid  computer/Unimoment-Monte  Carlo  area.  During  the  period  15  Sept. 

to  1 Nov. , 1978  a search  for  alternative  computational  methods  was 

undertaken  with  the  goal  of  finding  a method  which  would  allow  the  goals 

of  this  study  to  be  fulfilled  while  overcoming  the  difficulties  listed 

above.  Three  techniques  were  considered: 

1.  The  method  of  moments  [Harrington,  1968],  which  would  be  used 

to  calculate  electric  currents  induced  on  the  surface  of  the 

patch,  from  which  the  far  fields  could  be  found  using  the 

magnetic  vector  potential  method. 

2.  The  finite  difference  method,  in  which  differential  equations 
are  approximated  by  difference  equations.  These  are  solved  for 
the  electric  field  at  the  perimeter  of  the  patch,  and  the 

far  field  may  be  obtained  by  using  the  electric  vector  potential. 

3.  The  finite  element  method  [Strang  & Fix,  1973],  a variational 
method  in  which  the  minimization  process  automatically  seeks 
out  the  solution  which  is  "closest"  to  the  true  analytical 
solution.  The  far  field  is  calculated  by  the  same  method  as 
in  (2). 

Of  these  three,  the  finite  element  method  was  chosen  mainly  because  of 
the  generality  of  the  technique  (many  different  patch  shapes  may  be 
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considered)  and  the  mathematical  integrity  inherent  in  the  method 
(unlike  the  method  of  moments,  which  only  may  converge  - and  sometimes 
to  the  wrong  solution) . 

3. 1 Mathematical  Development  and  Numerical  Implementation 

Maxwell's  two  curl  equations  may  be  combined  to  form  the  Helmholtz 
equation  (with  zero  charge  density) : 


V2E  + k2  E = juuJ 


(3-1) 


where  k^  = 2tt/Xq  = u)/y  e ^ree  sPace-  The  Helmholtz  equation  itself 

does  not  couple  together  the  rectangular  components  of  E,  i.e. 


V2E  + kZ  E = iooyJ 

X o X X 

V2E  + k2  E = jioyJ 
y o y y 

V2E  + k2  E = jwyJ 
Z o z z 


(3-2) 


However,  the  boundary  conditions  may  couple  the  rectangular  components. 

For  example  the  radiation  condition  requires  that  the  radially  directed 
component  of  E (Er)  vanish  properly  at  infinity.  Since  Ef  is  composed 
of  (sin9  cos<t)Ex  + (sin0  sirnfOE^  + cos9Ez,  the  components  of  E are  coupled 
through  the  radiation  condition.  If  a vector  problem  is  to  be  avoided, 
then  the  boundary  conditions  specified  in  the  problem  must  contain  only 
one  rectangular  component  of  E.  The  report  by  Lo,  et.  al.  [1977]  assumed 
that  9Ez/3n  = 0 (open  circuit  at  radiating  edges,  or  Neumann  boundary 
value  problem)  and  obtained  fairly  good  results  for  a variety  of  patch 
shapes.  This  present  work  assumes  an  impedance  boundary  condition  on  E^ 
(Robbin  boundary  value  problem) : 


+ a E “0 
z 


(3-3) 


Finally,  it  is  assumed  that  Ez  does  not  vary  with  z.  This  is  a 
fairly  safe  assumption  for  the  lower-order  modes,  but  is  risky  for 
higher-order  modes.  By  making  this  assumption,  it  is  possible  to  reduce 
the  original  three-dimensional  problem  to  a two-dimensional  problem 
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of  a partial  differential  equation  (PDE)  plus  boundary  conditions  (BC) : 


+ 


32u 

"3 y 7 


2 

+ k“u 


jwpJz 


(PDE) 


3u 

3n 


+ a u 


- 0 


(BC  on  patch  perimeter) 


(3-4) 


where  for  simplicity,  the  variable  u has  been  substituted  for  E . 

z 

The  equivalent  variational  formulation  of  (3-4)  is  to  minimize 
the  functional  [Strang  & Fix,  pp.  70-71,  1973]: 


I(v) 


,2  2 

(v  + v 
x y 


interior 


k^v^  + 2jwuJzv)  dA  + 0 a v“  dZ 

perimeter 


(3-5) 


for  all  permissible  functions  v(x,y).  The  particular  function  v (x,y) 
which  minimizes  (3-5)  is  the  "best"  solution  to  (3-4).  Details  of  this 
development  may  be  found  in  Strang  & Fix  [1973]. 

Three  tasks  remain  in  order  to  solve  the  problem.  First,  a set  of 
"permissible"  v functions  must  be  determined.  The  goal  is  to  choose  trial 
functions  i(p  , . • • • <1^ , and  from  among  their  linear  combinations 

to  find  the  one  which  minimizes  (3-5).  In  theory,  there  always 
exists  a set  of  trial  functions  which  is  complete,  but  it  may  be  difficult 
or  even  impossible  to  compute  with  them.  The  underlying  idea  of  the  finite 
element  method  is  this:  divide  the  patch  into  smaller  pieces  such  as 
triangles  or  rectangles.  Then  within  each  piece  the  trial  functions  are 
given  a simple  form:  all  trial  functions  are  zero  except  one  which  is 
usually  a polynomial.  Boundary  (continuity)  conditions  are  imposed 
locally,  along  the  edge  of  a triangle  or  rectangle  rather  than  globally 
as  in  the  classical  Ritz  method.  The  approximation  of  the  solution  may  be 
refined,  if  necessary,  not  by  including  more  complex  trial  functions,  but 
by  further  subdivision  of  the  patch.  A digital  computer  is  ideally  suited 
to  this  task. 

Second,  the  integrals  in  (3-5)  must  be  computed  for  each  trial  function 
in  terms  of  still-to-be-determined  parameters.  For  the  case  at  hand  these 
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parameters  are  the  nodes  of  the  subdivisions,  but  for  other  physical 
problems,  the  parameters  may  represent  other  quantities.  The  integrals 
in  (3-5)  are  particularly  easy  to  compute  since  each  Vy  is  zero  except 
within  one  subdivision. 

Finally,  the  problem  must  be  put  into  a form  that  can  be  solved 
by  a computer,  as  hand  calculation  is  too  tedious  and  time-consuming. 
This  form  is  given  symbolically  by  [Strang  & Fix,  1973,  p.  32]: 


T T 

02  - 2 £ F 


(3-6) 


where  Q is  the  column  matrix  of  coefficients  (the  solution),  K is  the 
finite  element  matrix  arising  from  (3-5),  and  F is  the  forcing  function 
column  matrix,  the  linear  (in  v)  term  of  (3-5).  Equation  (3-6)  has  a 
minimum  for 


K£  = F 

and  (3-7)  is  the  matrix  problem  to  be  solved. 

The  eigenvalue  problem  (with  the  eigenvalue  A 


(3-7) 


-k  ) may  be  set  up  as 


K 2 = A K 2 


(3-8) 


where  is  obtained  from 


(v2  + v2)  dA  + 
x y 


interior 


a v dZ 


(3-9) 


perimeter 


and  is  obtained  from 


v2  dA 


interior 


(3-10) 


Notice  that  the  impedance  boundary  value  a does  play  a role  in  the  eigen- 
value determination  while  the  forcing  term  ju>uJz  does  not.  The  eigenvalues 
and  eigenvectors  of  (3-8)  may  be  obtained  by  using  standard  subroutine 
libraries  (SSP  or  IMSL  for  example  - see  list  of  references). 

The  calculation  of  the  and  K?  matrices  for  a general  polygonal 
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microstrip  antenna  has  been  implemented  in  the  computer  program  MICRO, 
a listing  of  which  can  be  obtained  from  the  authors  . The  important 
inputs  to  this  program  are  (1)  the  number  of  triangular  subdivisions 
made  to  the  original  patch,  (2)  the  numbers  of  the  three  nodes  which 
define  each  triangle,  (3)  the  total  number  of  nodes,  and  (4)  the 
x-y  coordinates  of  each  node.  From  this  data,  the  program  calculates 
the  and  matrices  for  the  Neumann  problem.  If  certain  nodes  are 
constrained  to  a field  of  zero  (Dirichlet  condition)  these  node  numbers 
are  also  entered.  The  eigenvalues  of  the  resulting  generalized  matrix 
problem  are  computed  and  checked  by  IMSL  subroutine  EIGZF  and  printed. 

3. 2 Numerical  Difficulties  Overcome  or  to  be  Overcome 
3.2.1  Matrix  Calculation  for  Irregular  Polygon  Patches 


To  utilize  the  full  power  of  the  finite  element  method,  the  computer 
programs  written  to  implement  the  technique  must  be  general  enough  to  handle 
a variety  of  geometries.  For  regular  geometries  such  as  the  rectangle, 
circle,  etc. , it  is  simple  enough  to  devise  computer  programs  that  will 
automatically  partition  the  regular  region  into  appropriate  subdivisions. 

On  the  other  hand,  for  general  polygonal  microstrip  antennas,  such  a program 
would  be  unwieldy  at  best  and  not  general  enough  at  worst.  This  is 
particularly  important  in  the  case  of  patch  antennas  because  it  is  the 
unusual  geometry  that  is  often  so  interesting.  For  example,  there  is  a 
five-sided  polygon  patch  which  can  be  excited  at  a single  feed  point  to 
give  right  circular  polarization,  at  a second  point  to  give  left  circular 
polarization,  and  at  a third  point  to  given  linear  polarization  [Weinschel, 
1975].  To  make  the  MICRO  program  as  general  as  possible  the  so-called 
area  method  [Strang  & Fix,  1973,  pp.  90-97]  was  used  rather  than  the 
node  method  more  familiar  to  those  who  work  with  moment  methods. 

The  area  method  works  as  follows.  Each  triangular  subdivision  area  is 

taken,  one  at  a time,  and  the  three  nodes  which  define  the  triangle  are 

determined.  The  integrals  in  (3-5)  are  performed  with  these  three  unknown 

node  values  as  parameters.  The  resultant  expression  is  a quadratic  and 

may  be  put  into  the  form  of  (3-6)  to  determine  what  is  called  K^,  the 

elemental  K matrix.  K.  is  different  for  each  area  and  contains  at  most 
— — e 

nine  non-zero  entries.  The  sum  of  all  K matrices  is  the  K matrix  of  (3-6). 

— e — 
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As  an  example,  suppose  one  triangular  area  is  bounded  by  nodes  1,  2,  and 
4.  Then  K has  the  form: 


X 

X 

0 

X 

X 

X 

0 

X 

0 

0 

0 

0 

X 

X 

0 

X 

(3-11) 


k 

i 
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In  this  matrix,  only  k. . / 0 for  i,j  = 1,2,  or  4.  A detailed  comparison 
between  the  area  method  and  the  nodal  method  is  given  in  Strang  & Fix 
[1973]. 

The  second  matrix  computational  difficulty  was  the  evaluation  of  the 
integrals  of  (3-5).  These  integrals  are  of  the  form 


triangle 


s 

y 


dA 


(3-12) 


Ho land  and  Bell  [1969]  have  given  explicit  formulas  for  these  integrals, 
keeping  the  origin  of  the  local  (x,y)  system  at  the  centroid  of  the  triangle. 
For  r + s < 6, 

Prs  ' cr«  A (X[Y?  + X[Y2  + X5Y?>  <3-13> 

where  (X^,Y^),  (X^jY^)  are  the  vertices  of  the  triangle  in  the 

centroid  system,  A is  the  area,  and  c^  = 1,  c^  = 0,  C2  = 1/12,  c^  = = 

1/30,  Cj.  = 2/105.  Other  possible  coordinate  systems  and  formulas  are  found 
in  Strang  and  Fix  [1973.,  pp.  90-96]. 

3.2.2  Impedance  Boundary  Conditions 

The  impedance  boundary  condition  is  easy  to  handle  numerically  once 
the  impedance  parameter  a(x,y)  of  equation  (3-4)  has  been  determined 
theoretically.  This  is  related  to  the  wall  admittance  per  unit  length. 
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Computationally,  terms  are  added  to  the  K matrix  corresponding  to  that 
portion  of  (3-5)  containing  o.  If  the  patch  antenna  has  been  divided 
into  sufficiently  many  sections  then  a may  be  assumed  constant  over  any 
one  section,  and  the  integration  is  trivial.  Otherwise  the  integration 
may  be  handled  numerically. 

One  significant  problem  must  be  dealt  with  whenever  a is  imaginary  or 
complex  (indicating  that  energy  is  being  radiated  from  the  cavity  walls) : 
the  K or  K^matrices  become  complex.  For  the  radiation  problem  this  presents 
no  difficulty  since  complex  arithmetic  may  be  used  in  the  matrix  inversion 
routine.  However,  for  the  eigenvalue/oigenvector  problem,  only  real  arith- 
metic routines  have  been  found.  Hence,  a conversion  from  an  N x N complex 
problem  to  a 2N  x 2N  real  problem  must  be  performed.  Mathematically,  this 
is  trivial  but  such  an  increase  in  matrix  size  may  require  over  times  as 
much  computer  time. 

3.2.3  Eigenvalue/ Eigenvector  Computation 

Up  to  this  point  all  eigenvalues  and/or  eigenvector  calculations  have 
been  made  by  using  IMSL  subroutine  EIGZF.  However,  to  increase  the  size 
of  the  problems  now  being  solved  or  to  include  complex  matrices  in  the 
subsequent  work  under  this  grant,  a different  approach  must  be  taken.  The 
most  promising  technique  thus  far  found  is  called  "subspace  iteration." 

It  is  possible  to  compute  the  first  p eigenvalues/vectors  in  about  8 
iterations  if  the  subspace  dimension  Z.  was 

Z = min(2p , p+8)  (3-14) 


Since  only  the  lower  few  modes  are  of  importance  in  microstrip  antenna 

work  (and  the  higher  modes  are  just  numerical  garbage)  this  technique  seems 

to  be  the  most  promising.  Briefly,  to  solve  A X = X M X , start  with  an 

eigenvector  guess  of  X .,  which  is  a matrix  with  Z orthonormal  columns. 

— n- 1 

Solve  K Y = M X , . Then  solve  the  eigenvalue  problem 
n n—  1 


(YT  K Y ) Q 
— n n 


v (Y  M Y ) 0 
—n  — — n 


(3-15) 


The  v,  are  the  approximate  eigenvalues  and  the  new  matrix  X of  approximate 
1 — n 

eigenvectors  is  formed  by  multiplying  Y^  with  the  ZxZ  matrix  of  eigenvectors 
of  (3-15).  Since  equation  (3-15)  is  of  only  dimension  Z,  subroutine  EIGZF 
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may  be  used  efficiently  here.  For  example,  suppose  the  first  four  modes 
are  needed  in  a microstrip  problem.  Then  by  (3-14),  t = 12,  and.  only 
about  eight  12  x 12  eigenvalue  problems  need  be  solved. 

3.2.4  Size  of  Matrices  versus  Core  Memory 

At  PSL  only  about  140K  bytes  of  computer  storage  are  avilable  for 
matrices  and  K^.  This  limits  the  size  of  and  K9  to  about  132  x 132 
four-byte  words.  However,  most  of  the  entries  in  either  matrix  are  zero. 
For  example,  in  a regular  triangular  partition  there  are  at  most  seven 
entries  per  row  regardless  of  the  size  of  the  problem  to  be  solved.  To 
utilize  more  effectively  the  available  storage  area,  a linked  list  sparse 
matrix  routine  has  been  written  following  the  guidelines  of  Wexler  [1978] 
and  Zollenkopf  and  Reid  [1971].  This  routine  has  been  tested  for  the 
simple  case  of  the  static  Dirichlet  problem  within  the  rectangle.  The 
test,  while  verifying  the  integrity  of  the  program,  also  illustrates  the 
tremendous  savings  in  computer  memory  that  can  be  obtained.  For  the 
test  case,  a ten-to-one  reduction  in  memory  was  achieved. 

One  of  the  advantages  of  using  a sparse  matrix  program  is  the  ease 
with  which  the  program  may  be  combined  with  subspace  iteration.  Subspace 
iteration  requires  only  txZ  matrices  (see  3.2.3)  where  t is  small.  Thus 
the  and  Kj  data  may  be  kept  in  a linked  list  rather  than  in  a full  dense 
matrix,  and  the  storage  requirements  may  be  kept  to  on  the  order  of  LyJL. 

As  an  example,  consider  the  finite  element  solution  of  the  rectangular 
interior  problem  for  which  one  hundred  (10  x 10)  interior  nodes  have  been 
specified.  The  lowest  three  modes  are  desired.  Using  traditional  matrix 
techniques,  two  100  x 100  matrices  would  have  been  needed  to  solve  for  the 
eigenvectors  and  eigenvalues  of  the  problem,  or  about  80K  bytes  of  memory. 
However,  all  but  460  entries  in  either  matrix  are  zero.  Sparse  matrix 
techniques  utilize  algorithms  in  which  none  of  the  zero  entries  need  be 
stored.  Instead,  two  auxiliary  tables  are  used.  For  the  above  example, 
less  than  7.6  K bytes  of  matrix  space  plus  3.8  K bytes  of  table  space 
would  be  required,  a savings  of  86%.  (In  addition,  solution  times  are 
generally  more  rapid,  offering  further  savings  in  computer  resources.) 

The  use  of  subspace  iteration  gives  a great  reduction  in  computer  time 
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because  only  the  lowest  three  modes  are  needed  in  this  example.  Using 

equation  (3-14) , & = 11.  Hence,  eight  11  x 11  eigen-problems  need  be 

solved  versus  one  100  x 100  eigen-problem.  Since  solution  times  vary  at 

2 

least  on  the  order  of  N , this  represents  a savings  of  90%.  On  the 
PSL  IBM  computer  this  represents  a savings  of  14  minutes  for  a 51  x 51 
problem. 

Finally,  it  should  be  noted  that  the  method  of  moments  generates 
in  general  full,  dense  matrices,  and  sparse  matrix  techniques  cannot  be 
used.  This  is  because  most  moment  method  problems  arise  from  integral 
equations  whereas  most  finite  element  method  problems  arise  from  partial 
differential  equations. 

At  this  writing  neither  the  sparse  matrix  routine  nor  the  subspace 
iteration  program  have  been  implemented  per  se  into  the  MICRO  program. 

3.3.  Handling  the  Feedpoint 

A microstrip  antenna  is  usually  fed  in  one  of  two  ways:  (1)  coaxially 
from  beneath  the  ground  plane,  or  (2)  on  the  edge  of  the  patch  by  a micro- 
strip transmission  line.  Each  method  of  feeding  may  be  handled  within  the 
framework  of  the  present  mathematics. 

The  coaxial  feed  may  be  modeled  as  a point  current  source  at  (x^y^) 
[vide  Sec.  2.3.5]: 

J (x,y)  = J <5(x  - x ) 6 (y  - y ) (3-16) 

z o o o 

and  may  be  treated  as  a forcing  term  in  equations  (3-4)  and  (3-5).  The 
integrals  for  the  F-matrix  of  (3-7)  are  simplified  because  of  the  Dirac 
delta  functions.  It  has  been  assumed  that  has  no  z-dependence  due  to 
the  relatively  small  spacing  (much  less  than  a wavelength)  between  patch 
and  ground  plane  on  most  microstrip  antennas.  This  assumption  may  not  be 
valid  for  feeds  placed  near  an  edge  or  horizontal  metallic  struction.  In 
these  cases,  a more  rigorous  feedpoint  attachment  theory  must  be  used 
[Newman,  1978]. 

The  microstrip  line  feed  may  be  handled  via  the  boundary  conditions. 
Over  the  width  of  the  line,  assume  that 

E (x,y)  *»  E (3-17) 

z o 
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and  elsewhere  use  the  impedance  boundary  condition  discussed  in  Sec.  3.1. 
The  effect  of  (3-17)  is  to  constrain  one  or  more  unknowns  to  the  value 
Eq;  this  has  been  discussed  previously  in  Sec.  2.3.5.  These  node  values 
are  then  removed  from  the  K matrices  and  transferred  to  the  _F  matrix 
as  equivalent  constants  (knowns) . While  this  formulation  requires  that 
the  boundary  conditions  be  mixed  (a  rather  intractable  theoretical 
problem),  mixed  boundary  conditions  are  not  at  all  difficult  to  handle 
numerically. 
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4.0  NUMERICAL  RESULTS 

4 . 1 The  Rectangular  Patch 

4.1.1  Closed  Cavity  Formulation  (Neumann  Boundary  Condition) 

Consider  a rectangular  microstrip  cavity  of  length  a and  width  b, 

as  shown  in  Fig.  4.1.  If  there  is  no  radiation  out  the  side  walls,  the 
component  of  the  magnetic  field  will  be  zero  on  the  y = 0 and  y = b 
walls;  the  H components  of  the  magnetic  field  will  be  zero  on  the  x = 0 
and  x = a walls.  The  mode  vectors  (c.f.  Sec.  2.3.5)  are  TM  to  z and 
for  this  Neumann  boundary  condition  are  given  by 
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2 if  m 0 and  n ^ 0 
Jl  if  m = 0 or  n 5 0 


(4-2) 


(4-3) 


which  satisfies  the  orthonormality  condition  (2-29).  The  corresponding 


resonant  wavenumbers  k are  given  by 

mn  ■ 


k = /(mir/b)2  + (nir/a) 2 
mn 

from  which  the  resonant  frequencies  may  be  obtained  by 

k 
f 


(4-4) 
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(4-5) 


The  magnetic  current  distribution  at  the  four  edges  is  computed  from 
(c.f.  eqn.  2-46) : 

M = z E x n 
s z 

where  E,  is  the  z-directed  electric  field  at  the  radiating  wall. 
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The  resonant  wavenumber  k may  also  be  computed  numerically,  using 
the  finite  element  method  outlined  in  Section  3 of  this  report.  To 
illustrate  the  comparison,  consider  a rectangular  microstrip  cavity  with 
a = 10  cm  and  b = 20  cm;  the  interior  dielectric  is  assumed  to  be  air. 
Table  4-1  compares  the  resonant  frequencies  calculated  using  (4-5)  to 
that  computed  from  the  method  of  finite  elements  using  100  square 
subdivisions.  Two  points  are  noted: 

1.  The  resonant  frequencies  computed  by  finite  elements  are 
greater  than  the  true  values. 

2.  As  the  mode  number  becomes  higher,  the  error  becomes  greater. 


Table  4-1 


Comparison  of  Computed  and  Theoretical  Resonant 
Frequencies  for  Rectangular  Microstrip  Cavity 
with  b = 2a  = 20  cm. 


Mode  Number 

(m,n) 

Resonant 

True 

Frequency 

Finite  Elements  Approx. 

(0,0) 

0 

0 

(1,0) 

750.00  MHz 

750.48  MHz 

(0,1) 

1500.00 

1502.26 

(2,0) 

1500.00 

1503.76 

(1,1) 

1677.05 

1720.83 

(2,1) 

2121.32 

2211.62 

(0,2) 

3000.00 

3268.77 

1 

Both  phenomena  are  characteristic  of  the  finite  elements  method.  However, 
for  the  microstrip  antenna  problem,  neither  is  of  concern  because  (1) 
sufficient  subdivisions  can  be  taken  to  insure  that  the  computed  resonant 
frequency  is  precise  enough  [Strang  & Fix,  1973,  Ch.  2]  and  (2) 
only  the  first  three  or  four  modes  are  of  practical  interest. 

Since  there  is  no  radiation,  the  mode  vectors  (4-1)  form  a complete 
orthonormal  set  of  discrete  eigenfunctions  (c.f.  Sec.  2.3.5)  which  can  be 
used  to  represent  any  arbitrary  interior  electric  field,  i.e. 
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This  double  Fourier  series  has  coefficients  which  may  be  found  using 
the  orthonormality  property  of  the  modal  vectors;  since  this  is  a 
straightforward  exercise  in  mathematical  physics,  only  the  result  will 
be  given.  Assuming  the  cavity  is  fed  by  a z-directed  filamentary  current 
located  at  (xQ,yo),  the  total  interior  electric  field  is  given  by 
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(4-8) 


If  the  cavity  thickness  is  small,  I(z')  =1.  ; furthermore,  the  input 

i ln 

voltage  is  = -t  E (xQ,yQ)  so  that  the  input  impedance  may  be  calculated 
from  (2-28) : 
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Thus,  for  the  rectangular  cavity 
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(4-10) 


The  input  impedance  is  purely  reactive  since  no  allowance  for  dissipated  or 
real  radiated  power  has  yet  been  made.  This  will  be  considered  in  the  next 
section. 

It  should  be  remarked  that  the  shape  and  polarization  of  the  far-field 
pattern  can  be  fairly  accurately  predicted  using  only  the  magnetic  line 
source  model,  i.e.  eqns.  (2-58)  through  (2-61),  with  magnetic  surface  currents 
computed  from  (4-6)  and  (4-1).  This  simple  approach  works  well  for  predicting 
the  pattern  shape  for  discrete  modes,  as  has  been  pointed  out  by  Lo,  et.  al. 
[1977]. 
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4.1.2  Radiating  Cavity  Formulation  (Impedance  Boundary  Condition) 

Since  the  microstrip  antenna  cavity  radiates  real  and  reactive 
power,  the  interior  fields  must  satisfy  an  impedance- type  boundary 
condition  at  the  radiating  walls.  If  the  wall  admittance  is  purely 
susceptive,  then  the  eigenvalues  k^  (similar  to  mr/a  for  the  discrete 
case)  and  k (similar  to  mrr/b  ) will  still  be  real-valued  and  discrete. 
However,  for  the  radiative  condition,  the  wall  admittance  is  complex. 

It  may  be  shown  that  the  eigenvalues  are  given  by  the  coupled  trans- 
cendental equations: 
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(4-11) 


(4-12) 


(4-13) 


(4-14) 

(4-15) 


with  Y being  the  wall  admittance  taken  on  the  appropriate  wall;  Z = 376  ft 
w o 

and  Xq  is  the  free-space  wavelength. 

In  the  dominant  mode,  effectively  only  the  y = 0 and  y = b 

walls  radiate;  furthermore,  the  electric  field  in  the  radiating  walls  is 
essentially  independent  of  x so  that  the  wall  conductance  may  be  approximated 
by  that  of  a parallel-plate  TEM  waveguide  opening  onto  a large  ground  plane; 
the  resulting  conductance  is  given  by  [Harrington,  1961,  p.  183]: 


Gw  " 376  fl  - ^5(t/Xo^[ano]  (V) 

For  electrically  small  substrate  thicknesses,  this  becomes 


(4-16) 


1 
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(4-17) 


G = 0.00836  (a/X  ) (U) 

w o 

This  is  plotted  in  Fig.  4.2  along  with  the  conductance  of  a slot  antenna 
radiating  into  half-space  (Derneryd,  1976);  it  is  seen  that  the  conductance 
computed  from  (4-17)  is  roughly  double  that  used  by  Derneryd  for  a slot 
antenna. 

The  wall  susceptance  may  be  computed  from  an  expression  for  the 
capacitance  of  an  open  microstrip  circuit  [Hammerstad,  1975;  Derneryd,  1976] 
and  may  be  written  in  the  convenient  form 

kJL  a 

B = 0.01668  — — e OJ)  (4-18) 

w t A e 

o 


where 


a 

— = 0.412 

e + 
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.300 

r—  + 
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.258 

a 

.813 

t 

e 

— + 
l t 

(4-19) 


and  where  is  the  effective  dielectric  constant  given  by  [Schneider,  1972]: 


e = 
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e +1 
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e - 1 
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lOt 
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(4-20) 


Fig.  4.3  is  a graph  of  both  the  susceptance  computed  from  (4-18)  and  that  of 
the  parallel-plate  waveguide  radiating  into  a half-space;  both  are  for  an 
assumed  air  dielectric. 

For  the  TM-^qq  mode,  = 0 and  the  resonant  frequency  may  be  computed 


from 


k 


f = 4.775  — (GHz)  (4-21) 

r JT 

e 

where  k (in  cm  'S  is  obtained  from  a solution  to  the  transcendental  equation 

y 

(4-12)  and  where  a is  computed  from  (4-15).  The  effect  of  the  impedance 
boundary  condition  is  to  produce  a dependence  of  the  resonant  frequency  on 
the  dimension  a;  this  is  not  found  in  the  Neumann  condition  (c.f.  eqns.  4-4 
and  4-5  with  n = 0) . This  theory  produces  good  agreement  with  experimentally- 
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SUSCEPTANCE 
( MHOS’) 


MICROSTRIP  EDGE 
EFFECT-CAPACITANCE 


AIR  DIELECTRIC  ASSUMED 


a/' A = .50 


PARALLEL- PLATE 
WAVEGUIDE  END  EFFECT 


PARALLEL-PLATE 


a/\  = .50 


MICROSTRIP 


SLOT  WIDTH  / WAVELENGTH  t/A 


Fig.  4.3  Susceptance  as  a function  of  slot  width  for  a microstrip 
open  circuit  and  for  a parallel-plate  waveguide  radiating 
into  half-space.  An  air  dielectric  is  assumed  for  both 


cases. 
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measured  resonant  frequencies  at  both  L-band  [Lo,  et.al.,  Wy]  and  S-band, 
as  is  shown  in  Figs.  4.4  and  4.5.  In  both  cases,  it  is  noted  that  as  the 
dimension  a is  increased  the  resonant  frequency  decreases;  this  can  be  seen 
more  clearly  in  Fig.  4.6  which  is  a graph  of  both  left-  and  right-hand 
sides  of  (4-12).  The  solution  for  the  TM^qq  mode  in  the  Neumann  condition 
is  k^b  = it;  however,  with  an  impedance  condition,  this  becomes  slightly  less, 
typically  96%  to  98%  of  it. 

The  values  of  a are  complex,  although  the  real  part  is  normally 
much  larger  than  the  imaginary  part  (because  of  the  susceptance  being  larger 
than  the  conductance).  Thus,  in  solving  (4-15),  only  the  real  part  of  a is 
used  and  the  eigenvalues  are  approximately  real. 

The  interior  fields  and  input  impedance  may  be  calculated  using  a 
modification  of  the  Fourier  series  expansion  presented  in  Sec.  4.1.1  by 
letting  the  resonant  frequencies  become  complex  [Harrington,  1961,  p.  435]: 


“2r  <!  + 


(4-22) 


where  is  the  TM  resonant  frequency.  In  the  vicinity  of  resonance, 
r mno  ^ 1 

(4-10)  becomes 


miry  nirx 

2 ? O 0 O 

wt  Y_  cos'2  — — cos^  

mn  b a 


(D2  - co2  (1  + *!-  ) 
mn  Q 


(4-23) 


where  is  the  series  inductance  of  the  probe  feed  connected  to  the  micro- 
strip patch,  i.e. 


( n ) 


(4-24) 


The  second  term  of  (4-23)  is  of  the  same  form  as  the  input  impedance  to  a 
parallel  RLC  circuit  as  shown  in  Fig.  4.7;  thus, 


Q A t miry  nirx 

mn  i o o o 

R r—  cos'2  — cos'2  

uj  cab  b a 

mn 


1 y ^ t miry  nirx 

mn  ■>  o o o 

; — COS*  — COS*  

2 cab  b a 


(4-25) 


(4-26) 


* It  is  assumed  here  that  only  one  resonant  mode  is  excited,  i.e.  no 
closely-spaced  modes  or  degeneracies. 
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Fig.  4.5  A comparison  between  experimentally-determined 
resonant  frequencies  and  theoretical  values 
based  on  an  impedance  boundary  condition. 

S-band  microstrip  patch  fed  coaxially;  experimental 
values  include  compensation  for  probe  series 
inductance  effect. 
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Lp  = SERIES  INDUCTANCE  ASSOCIATED  WITH  PROBE  FEED 

R = RESISTANCE  ASSOCIATED  WITH  RADIATION  AND 
HEAT  LOSSES 

L = INDUCTANCE  ASSOCIATED  WITH  ENERGY  STORAGE  IN 
CAVITY 

C = CAPACITANCE  ASSOCIATED  WITH  ENERGY  STORAGE  IN 
CAVITY 


Fig.  4.7  Equivalent  circuit  model  of  microstrip  cavity 
fed  by  coaxial  probe. 
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(4-27) 


Also  it  is  noted  that 


(4-28) 


u)  R C 
run 


(4-29) 


The  resistance  may  be  computed  by  noting  that  at  the  resonant  frequency 
(where  the  second  term  in  4-23  becomes  purely  real),  the  input  resistance  is 


R.  = R + R,  + R 

in  c d 


(4-30) 


where  Rc  is  the  equivalent  resistance  associated  with  dissipative  losses  in 
the  conductors  (c.f.  eqn.  2-17)  and  R^  Is  the  equivalent  resistance  associated 
with  dissipative  losses  in  the  dielectric;  R is  then  the  reciprocal  of  twice 
the  conductance  of  one  slot  (eqn.  4-17),  i.e.  for  the  dominant  mode 


R = 59.8  (A  / a)  ( Q ) 

o 


(4-31) 


The  first  term  in  (4-27)  will  be  recognized  as  the  d.c.  capacitance  of  a 
parallel-plate  capacitor  of  area  ab  and  thickness  t. 

It  is  emphasized  that  in  (4-23)  the  series  inductance  term  is 
associated  with  the  coaxial  probe  and  results  from  coupling  to  the  TM  field. 
For  a microstrip  antenna  fed  by  a microstrip  feedline  this  term  is  not 
present,  and  the  input  impedance  is  given  by  the  second  term  of  (4-23). 

Fig.  4-8  compares  the  input  impedance  of  the  present  model  with  the 
transmission-line  model  of  Derneryd  [1976]  for  a microstrip-fed  patch  with 
a = 6.858  cm,  b = 4.140  cm,  t = 0.1588  cm  and  = 2.62.  Two  significant 
differences  are  noted.  First,  the  input  resistance  at  resonance  from  (4-23) 
is  118  Q whereas  that  from  the  Derneryd  model  is  220  0;  this  results  from 
the  different  models  used  for  the  slot  conductance.  Second,  the  resonant 


w*  T. 

«*> 


50 


I 

I 

I 


frequency  for  the  modal  expansion  model  is  predicted  to  be  2194  MHz, 
whereas  the  Derneryd  model  predicts  2172  MHz.  The  measured  resonant 
frequency  for  this  antenna  was  2195  MHz,  in  nearly  perfect  agreement  with 
the  modal  expansion  model. 

For  the  ™^qq  mode,  the  impedance  is  predicted  to  be  independent  of 
the  probe  position  xq;  this  is  verified  by  experimental  data  [Weinschel  & 

Lanphere,  1973]  for  feed  points  ranging  from  the  center  to  near  the  corner. 

When  the  patch  is  fed  by  a probe  (coaxially),  the  inductive  term  must 
be  retained  in  (4-23),  which  has  the  effect  of  moving  the  impedance  circle 
toward  the  inductive  side  of  the  Smith  Chart.  For  a typical  substrate  with 
er  = 2.62  and  t = 0.1588  cm,  = 7.7  !1  at  1 GHz,  15.5  0 at  2 GHz,  and 
80.1  ft  at  10  GHz.  Since  practical  impedance  measurements  are  often  made 
with  coaxially-fed  patches,  it  is  important  to  correct  such  measurements 
for  the  inductive  effect,  particularly  at  S-band  or  higher. 

As  the  feed-point  yQ  is  moved  inward  toward  the  center  of  the  element, 
the  resistance  decreases  in  accordance  with  (4-25) ; the  capacitance  increases 
but  the  inductance  decreases  in  such  a way  as  to  maintain  the  same  resonant 
frequency  and  the  same  Q. 

Figs.  4.9  - 4.12  compare  measured  impedance  data  with  theoretical 
values  (based  on  4-23)  for  a patch  with  a centered  feedpoint  at  several 
inset  distances  ranging  from  0 cm  (edge)  to  1.524  cm.  The  agreement  is  quite 
good,  and  seems  to  validate  many  of  the  preceding  assumptions.  Fig.  4.13 
compares  experimental  [Dunn,  1978]  and  theoretical  impedances  for  an  edge-fed 
L-band  patch  and  again  the  agreement  is  good. 

One  restriction  on  the  present  formulation,  which  assumes  a single 

I 

well-isolated  resonant  mode,  is  that  it  becomes  less  accurate  as  the 
dimension  a becomes  less  than  b.  Thus,  the  second  term  in  (4-23)  includes 
only  one  mode;  future  formulations  will  include  multiple  modes  and  provisions 
for  estimating  the  effect  of  feed  probe  slight  positioning  errors  on 
degenerate  modes  (e.g.  TMQ10  and  ™ioO  for  3 scluare  or  nearly-square  patch). 
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Fig.  4.12  Comparison  between  theoretical  and  experimental 
values  of  impedance  for  S-band  patch  fed  at 
1.524  cm  in  from  edge. 


The  power  dissipated  in  the  rectangular  patch  may  be  divided  into 

a portion  associated  with  conductor  losses  and  the  remainder  caused  by 

dielectric  losses,  as  described  in  eqns.  (2-16)  - (2-18).  By  using  the 

preceding  modal  expansions,  it  is  easy  to  show  that  at  a well-isolated 

resonant  frequency,  the  power  lost  as  heat  in  the  conducting  walls 

* 

(assumed  to  be  copper)  is  : 


where 


P.  = |l  |2R 
dc  1 in 1 c 


R = 0.00027/F”  — Q ( ft  ) (f  = frequency  in  GHz) 
c 3. 


(4-32) 


(4-33) 


The  power  lost  as  heat  in  the  dielectric  is 


P = |l  |2R 
dd  in  d 


where 


R , = 30  q2 

d e ab 

r 


( a ) 


(4-35) 


where  is  the  real  part  of  the  relative  permittivity  and  tan  5 is  the 
loss  tangent.  The  input  power  is 


I.  |2R. 
in  in 


(4-36) 


The  radiation  resistance  referred  to  the  input  is  therefore 


R = R.  - R - R , 
m c d 


(4-37) 


where 


R = 59.8  (1  /a)  ( 0,  ) 

o 


(4-38) 


Also,  Q - = ojW/ (Re  P^)  where  Re  = G(j|v|2  and  W = ^C^,|v|2  and  where 


CT  includes  both  the  patch  and  edge-effect  capacitance.  This  gives 


Q = (erb/2t)  + (2e£A£/t) 
Thus,  the  antenna  efficiency  is 


q = 1 - 


P + P 
dc  dd 


R + R , 
c d 


(4-39) 


(4-40) 


It  is  assumed  that  the  patch  is  fed  on  the  edge  and  in  the  THinn  mode. 
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^ As  a numerical  example,  let  a = 6.858  cm,  b = 4.140  cm,  t = .1588  cm. 


e'  = 2.62  and  X =13.6 
r o 

cm.  Then,  from  (4-33),  R =.38  Q and 

c 

from 

(4-38) , we  get  Q = 37, 

R =119  Q.  Table  4-2  shows  the 

effect 

of  various 

loss  tangents  on  both  R, 

(1 

and  the  antenna  efficiency  . 

Table  4-2 

Antenna  Efficiencies  for  Microstrip 

Patch 

Loss  Tangent 

Equiv.  Diel.  Resistance 

Efficiency  (Loss) 

tan  6 

Rd 

- n 

0.0001 

0.12  ft 

99.58 

% (.02  dB) 

0.0005 

0.60 

99.18 

(.04  ) 

0.001 

1.19 

98.68 

(.06  ) 

0.005 

5.96 

94.67 

(.24  ) 

0.01 

11.92 

89.66 

(.47  ) 

Fig.  4.13-1  is  a graph  of  the  antenna  efficiency  as  a function  of  substrate 
thickness  for  several  loss  tangents  and  for  the  microstrip  patch  in  the 
above  example. 
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Fig.  4.13-1  Antenna  efficiency  (gain/directivity)  vs.  substrate 
thickness  for  a rectangular  patch  S-band  microstrip 
antenna . 
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4.1.3  Multiple  Mode  Effects 


In  the  preceding  formulation  of  (4-23)  for  the  impedance  of  a 
rectangular  patch,  it  was  assumed  that  there  was  essentially  one  resonant 
mode  whose  bandwidth  was  much  less  than  the  frequency  separation  to  the 
next  adjacent  modes.  However,  for  some  microstrip  antenna  designs,  this  is 
not  necessarily  the  case;  for  example,  a square  patch  fed  at  a corner 
can  excite  both  the  TM^qq  and  TM^.^  modes  in  phase  and  at  ‘he  same 
resonant  frequency. 

Consider  a corner-fed  rectangular  patch  where  the  dimensions  a and  b 
are  approximately  the  same.  The  total  input  impedance  is  the  sum  of  the 
impedances  of  each  mode;  neglecting  the  inductive  reactance  of  the  probe  and 
using  only  the  second  term  of  (4-23)  with  xq  = yQ  = 0,  the  input  impedance  is 


Zin  Z100  + Z010 


where 


-2jwt 


'100 


eab  [w2  - (u2  (1  + ) ] 


10v  Q 


10 


-2jut 


'010 


eab  [w2  - w2  (1  + ^ ) 

ui  mq1 


(4-40a) 


(4-40b) 


(4-40c) 


It  is  assumed  that  all  other  modes  are  short-circuited. 

This  model  raises  the  possibility  of  increasing  the  bandwidth  of  the 
patch  by  stagger-tuning  the  modes.  The  resonant  frequency  is  determined 

primarily  by  the  dimension  b and  the  resonant  frequency  is  determined 

primarily  by  the  dimension  a. 

Figs.  4-14  a,b,  and  c show  this  effect  by  comparing  the  impedance 
behavior  of  square  and  nearly-square  patches,  all  with  a = 4.14  cm, 
t = 0.1588  cm  and  er  = 2.62.  Fig.  4-14a  shows  the  plot  for  b * a = 4.14  cm 
(square  patch),  which  is  the  degenerate  case.  Fig.  4-14b  is  for  b = 4.328  cm 
and  shows  a loop  which  is  associated  with  the  frequency  band  where 
and  are  about  the  same.  Fig.  4-14c  is  for  b = 4.500  cm  and  the  loop  is 

seen  to  become  larger  at  the  same  time  the  resonant  frequency  decreases. 

This  multiple-mode  excitation  technique  may  hold  promise  for  some  improvement 
in  bandwidth. 
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Multiple-mode  excitation  can  also  produce  various  polarization  states 
of  the  antenna.  For  a square  corner-fed  patch,  the  polarization  is  linear 
with  the  E-plane  oriented  along  the  diagonal  containing  the  feed.  For 
corner-fed  rectangular  patches  with  a t b,  however,  the  possibility  exists 
of  producing  circular  polarization  with  a single  feedpoint. 

By  using  (4-8)  and  (4-122),  the  interior  fields  of  a corner-fed 
rectangular  patch  may  be  written  as 


E(x,y)  = z 


JU)I 


in 


eab 


E E 
m,n 


miry  mrx 

cos  — cos  

b a 


j2 

mn 


(4-40d) 


mn 


For  a patch  with  dimension  a slightly  less  than  b,  the  and  ™Q10 

modes  will  have  closely-spaced  resonant  frequencies  so  that  only  two  terms 
of  the  series  are  retained: 
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(4-40e) 


The  magnetic  surface  current  on  the  y = 0 wall  is  thus  given  by 

TTX 
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and  the  magnetic  surface  current  on  the  x = 0 wall  is  given  by 
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(4-40f ) 


(4-40g) 


The  second  term  in  (4-40f)  and  the  first  term  in  (4-40g)  have  no  net  effect 
on  the  polarization  in  the  broadside  direction. 

To  achieve  circular  polarization,  lMgxi  = |M  | an<^  their  phases  must 
differ  by  90°. 

By  using  (4-40f)  and  (4-40g)  it  is  found  through  trial-and- error 
numerical  evaluation  that  a corner-fed  patch  with  b/a  = 1.029  should  produce 
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approximately  circular  polarization  in  the  broadside  direction.  A similar 
prediction  has  been  made  by  Kaloi  [1978]  although  no  ratio  of  b/a  was  given 
and  no  computed  or  measured  axial  ratios  were  presented.  The  numerical 
procedure  for  the  prediction  of  circular  polarization  is  to  enforce  the 
above  conditions,  i.e. 


and 


co 


“01 (1 


tan 


-1 


-M2n/Q 


10'  ^10 
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-1 


01/Q01 


= + (tt/2) 


00 


■ 00 


10 


01 


(4-40h) 


(4-40i) 


As  an  approximation  to  the  resonant  frequencies  oo^q  and  co^^,  the  Neumann 

boundary  condition  can  be  used  so  that  - Tr/b/ye”  and  = ir/a/ye. 

Then  = w1gR^oClO  where  R]_o  ~ 59.8AQ/a  and  = eab/2t;  similarly, 

Qnl  = to  .R  C where  R = 59. 8A  /b  and  CA1  = C1A.  These  equations  are 
01  01  01  01  01  o 01  10 

easily  solved  by  iterative  means  using  a programmable  calculator. 

As  an  example,  let  a = 4.14  cm  and  b = 4.26  cm.  The  predicted  on-axis 

axial  ratio  vs.  frequency  is  shown  in  Fig.  4-15.  If  the  patch  is  fed  at 

the  x^  = 0,  y^  = 0 corner,  the  polarization  is  left  circular;  if  the  patch 

is  fed  at  the  xq  = 0,  y^  = b corner,  the  polarization  is  right  circular. 

The  computed  polarization  bandwidth  (defined  as  the  frequency  spread  between 

axial  ratio  values  3 dB  greater  than  the  minimum)  is  22  MHz. 

A corner-fed  nearly-square  patch  with  the  above  dimension  was  constructed 

to  test  the  validity  of  the  above  theory.  The  dielectric  constant  was  2.62, 

the  thickness  was  1.588  mm  and  the  patch  was  fed  at  the  x = 0,  y =0  corner 

o o 

by  a coaxial  probe  (SMA  connector).  The  measured  axial  ratio  vs.  frequency 
is  shown  in  Fig.  4.15.  The  measured  impedance  diagram  is  shown  in  Fig.  4.15-1 
First,  the  agreement  between  theoretical  and  experimental  axial  ratios 
is  very  good.  Second,  it  is  observed  that  the  excitation  of  the  two  modes  is 
in  evidence  from  the  presence  of  the  cusp  (small  loop)  at  2200  MHz.  It  is 
also  noted  that  the  cusp  occurs  where  the  axial  ratio  is  minimum.  The  extensi< 
of  this  principle  relating  polarization  and  impedance  for  multi-mode  micro- 
strip patches  is  straightforward. 
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4 . 2 The  Circular  Patch 

Consider  a circular  microstrip  cavity  of  radius  a and  substrate 
thickness  t,  assumed  to  be  electrically  small.  The  orthonormal  electric 
field  mode  vectors,  which  are  again  assumed  TM  to  z , satisfy  the  Neumann 
condition  for  a non-radiating  cavity  and  are  given  by 


e.__(p,<t>)  = z J (k*  — ) cos  n<t> 

nr  a/2l7T  A n nr  3 

nr 


(4-41) 


where 


/2  if  n = 0 


if  n = 1 


(4-42) 


and 


J (k*  ) if  n = 0 
o or 


nr 


1 1/2 


I (n+l+2k)  in+1+2kOc;r) 
k=0 


(4-43) 


n i4  0 


t hi 

and  where  the  eigenvalue  k'  is  the  r -ranked  root  of  the  Neumann  condition 

nr 


J'  (k'  ) = 0 
n nr 


(4-44) 


Table  4-3  lists  in  ascending  order  of  frequency  the  eigenvalues  for  a 
few  of  the  lower-order  modes  [Abramowitz  and  Stegun,  1964]. 


Table  4-3 


Eigenvalues  of  J' (k'  ) = 0 
n nr 


Mode 

(n,  r) 

Eigenvalue 

k’ 

nr 

0,0 

0 

1,0 

1.84118 

2,0 

3.05424 

0,1 

3.83170 

2,1 

4.20119 

1,1 

5.33144 

2,1 

6.70613 

0,2 

7.01558 

r 

<r 


67 


The  resonant  wavelengths  are  then  calculated  by 


X = 

o 


2uV  e a 

r 

k' 

nr 


(4-45) 


Radiation  patterns  for  the  mode  have  been  measured  and  calculated 

by  Carver  [1976]  and  more  recently  and  completely  by  Lo,  et.  al.  [1977]. 
The  dominant  TM^qq  mode  has  a resonant  frequency  given  by 
c (1.84118) 


2ira 


(4-46) 
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where  c is  the  speed  of  light  in  vacuum.  The  normalization  constant  A 
for  this  mode  is  computed  numerically  from  (4-43)  to  be  0.33887.  The 
e^g  mode  is  zero  at  the  center  of  the  patch  and  has  a maximum  amplitude 
at  d>  = 0 and  i>  = it;  however,  the  electric  fields  in  the  edge  walls  are 
180°  out  of  phase  at  these  diametrically  opposite  points.  Such  a mode 
produces  a broadside  pattern  at  the  resonant  frequency  and  has  been  used 
successfully  in  practical  array  designs  [Parks  & Bailey,  1977].  This  mode 
may  be  excited  by  either  a probe  or  a microstrip  transmission  line  feed  at 
the  edge. 

As  before,  the  total  field  of  a probe-excited  cavity  can  be  expanded 
in  terms  of  the  mode  vectors.  For  a radiating  cavity  with  a finite  Q, 
this  becomes 


® 00  c e (p  ,<J>  ) e (p,<j>) 

E (p , (j>)  = - I (Jo>)  Z Z — ~ — ° 

z in  ,.,2  _ ,.,2  n 4-  1 


(4-47) 


n=0  r=0 


<r(1  + t > 


where  (P0’<t,0)  is  c^e  feed-point  location.  The  input  impedance  is  then 


00  00 


'in 


j to  E E 


C enr(po’<t>o) 


n-0  nmr  w2  - oj2  (1  + ^ ) 
nr  Q 


(4-48) 


In  the  vicinity  of  a resonant  frequency,  this  becomes 


'in 


J *L  * jW 


enr(po»V 

“2  - “^r(1  + Q> 


(4-49) 
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Using  A1 n = 0.33887  along  with  eqns.  (4-41)  and  (4-42),  the  mode  vector 
for  this  dominant  mode  is 

ein(p  .♦  ) - — j (1.84118  ^ ) cos  $ (4_50) 

The  first  term  in  (4-49)  is  the  series  inductive  reactance  associated  with 
the  probe  and  is  computed  from  eqn.  (4-24) . For  a circular  patch  fed  by 
a microstrip  transmission  line  this  first  term  is  not  present.  The  second 
term  is  of  the  same  form  as  the  impedance  of  an  RLC  circuit,  i.e. 
joj/C 

Z = — (4-51) 

“2-a,nr(1+Q) 


where 

aT  = 1/  (LC) 

nr 

Q = R/(u>L) 


(4-52) 

(4-53) 


Thus,  substituting  (4-50)  into  (4-49)  and  noting  the  analogies  (4-51)  - 
(4-53),  the  following  relationships  are  obtained: 
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Q 2. 775t 


nr 
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J2(l. 84118 
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(4-54) 
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1 2 . 775t  , p 

L 2-  jf(l. 84118  — ) (4-55) 

0)  ea  1 a 

nr 

ea2  p 

C = jT  (1.84118—)  (4-56) 

2.775t  1 

At  resonance,  the  input  resistance  is  equal  to  the  reciprocal  of  the  total 

effective  wall  conductance  G . An  exact  solution  to  this  theoretical 

w 

problem  of  the  wall  admittance  for  this  geometry  has  not  been  found  by  the 
author;  however,  an  approximate  form  will  be  given  in  a subsequent  paragraph. 

For  a radiating  cavity,  the  Neumann  condition  is  not  satisfied;  the 
impedance  boundary  condition  means  that  (1)  the  eigenvalues  and  eigenvectors 
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are  no  longer  purely  real,  and  (2)  the  resonant  frequencies  will  be  slightly 
less  than  those  obtained  for  the  closed  cavity  case.  The  transcendental 
equation  for  the  eigenvalue  of  the  mode  can  be  shown  to  be 

Ji(kio)  = - a Ji(kio)  <4~57> 

where 

Z Y t 
. o w 

a = j X i f4-58) 

o 

and  where  differentiation  is  with  respect  to  p,  Z = 376  fi,  and  Y is  the 

o w 

wall  admittance.  If  Y = 0,  then  (4-57)  reduces  to  the  Neumann  condition 

w 

(4-44)  for  n = 1.  By  use  of  a recursion  relationship,  (4-57)  can  be  written 
in  a more  convenient  form  as 

kio  vkio  > ■ <i  - “’w  <4-59) 

To  solve  this  equation  for  the  eigenvalue  , it  is  convenient  to  use 
the  polynomial  approximations  for  Jq  and  as  given  in  Abramowitz  and 
Stegun  [1961];  the  appropriate  value  will  be  slightly  less  than  the  Neumann 
condition  value,  i.e.  1.84118.  Table  4-4  lists  both  the  eigenvalues  and 
the  fractional  reduction  in  resonant  frequency  for  several  values  of  aa. 


Table  4-4 


Eigenvalues  for  Impedance  Boundary 
Condition  and  Fractional  Reduction 
in  Resonant  Frequency 


aa 

kio 

f / f 
r ro 

0 (Neumann  condition) 

1.84118 

1.00000 

-.01 

1.83345 

0.99580 

-.02 

1.82567 

0.99158 

-.02025 

1.82567 

0.99145 

-.03 

1.81780 

0.98730 

-.04 

1.80988 

0.98300 

-.05 

1.80190 

0.97867 
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The  relationship  of  the  values  from  (4-58)  to  the  values  tabulated  can 
be  made  directly  if  it  is  assumed  that  the  wall  susceptance  is  much  greater 
than  the  wall  conductance;  this  is  usually  the  case. 

Since  the  power  in  the  radiating  walls  varies  as  cos2!))  for  the  TM  qq 
mode,  the  effective  radiating  length  is  half  the  circumference.  Thus,  if 
a is  not  too  small,  the  relationships  for  a planar  parallel-plate  TEM 
guide  radiating  into  a half-space  (eqns.  4-17  and  4-18)  may  be  adapted  to 
the  circular  patch,  replacing  a (in  the  previous  equations)  with  Tra  in 
the  present  case,  so  that  the  wall  conductance  is 

Gw  * 0.00836  -p  (IT)  (4-60) 

o 

The  wall  susceptance  is  adapted  from  the  open  ended  microstrip  capacitor 
(eqns.  4-18)  and  is  thus 

Bw  = 0.00834  er  ( TJ  ) (4-61) 

o 

where  &£/ t has  been  taken  as  0.5. 

As  a numerical  example,  consider  a circular  patch  with  a = 6.75  cm, 

t * 0.1588  cm,  and  = 2.65.  The  resonant  frequency  for  the  Neumann 

condition  can  be  found  from  (4-46)  as  800.0  MHz.  From  (4-61),  B = .0125  IT 

-1  W 

and  from  (4-60),  G = 0.00178U.  From  (4-58),  a = -.003  cm  . A solution 
w 

to  the  impedance  boundary  condition,  with  aa  = -.02025,  gives 
k}  ; » 1.82567  as  shown  in  Table  4-4.  This  corresponds  to  a 99.145% 
reduction  in  the  resonant  frequency,  so  that  the  predicted  resonance  is 
793.2  MHz. 

This  example  was  chosen  to  correspond  to  measured  data  for  a patch 
of  the  same  size  [Lo,  et.  al.,  1977,  p.  36].  The  measured  resonant 
frequency  was  794  MHz,  which  is  in  nearly  perfect  agreement  with  the 
above  calculation.  This  appears  to  confirm  that  (4-61)  is  a good 
approximation.  Fig.  4.16  is  a comparison  between  the  theoretical 
input  impedance,  based  on  eqns.  4-49  (without  the  probe  inductance) 
through  4-61,  and  experimental  values  measured  by  Lo,  et.  al.  There  is 
generally  good  agreement,  except  that  the  measured  input  resistance  at 
resonance  is  slightly  larger  than  the  theoretical  value. 
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4 . 3 The  Pentagon  Patch 


In  Sec.  4.1.3,  it  was  pointed  out  that  circular  polarization  is 
theoretically  possible  from  a microstrip  antenna  with  a single  feedpoint, 
if  two  spatially  orthogonal  modes  were  excited  in  phase  quadrature. 

A practical  example  of  such  an  antenna  has  been  reported  by  Weinschel 
[1975]  who  obtained  circular  polarization  over  a 1%  bandwidth  from  a five- 
sided microstrip  patch,  shown  in  Fig.  4.17.  This  antenna  has  been  used  in 
an  S-band  cylindrical  array  producing  circular  polarization  [Weinschel,  1975] 
and  in  a UHF  scanning  array  for  buoy  or  ship  communication  to  a satellite 
[Weinschel  and  Carver,  1976]. 

A typical  impedance  plot  for  an  S-band  pentagon  microstrip  patch  is 
shown  in  Fig.  4.18.  The  loop  at  2279  MHz  is  evidence  of  the  existence  of 
two  modes,  as  discussed  previously  in  Sec.  4.1.3. 

The  modal  expansion  for  the  interior  fields  of  such  a polygon  cannot 
be  carried  out  by  the  classical  technique  of  separation  of  variables  since 
the  boundary  conditions  cannot  be  imposed  in  a simple  fashion.  The  method 
of  finite  elements  can  be  used,  however,  by  application  of  the  principles 
discussed  in  Sec.  3.0.  In  fact,  this  antenna  serves  as  an  illustration  of 
the  true  power  of  the  finite  element  method. 

Table  4-5  lists  the  computed  eigenvalues  (using  program  MICRO)  for 
the  first  eight  modes  of  a UHF  pentagon  patch  whose  first  non-zero  resonant 
frequency  was  measured  to  be  469  MHz.  The  grid  and  nodes  used  for 
computational  purposes  is  shown  in  Fig.  4.19. 

Table  4-5 

Results  of  Finite  Elements  Computed 
Values  of  Lower  Eight  Eigenvalues 
and  Resonant  Frequencies  for  Penta- 
gon Microstrip  Antenna 


n 

K 

n 

f 

n 

i 

0 -i 

0 

2 

0.1046  cm 

499  MHz 

3 

0.1112 

531 

4 

0.2902 

1384.4 

5 

0.4595 

2062.2 

6 

0.4799 

2153.6 

7 

0.7733 

3470.2 

8 

0.8672 

3891.3 
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The  resonant  frequencies  were  computed  assuming  a Neumann  boundary 
condition  (open  circuit  at  walls)  so  that  the  computed  values  are 
expected  to  be  a few  percent  high.  The  measured  f^  resonance  is  0.94 
times  the  calculated  value. 

Fig.  4.20  is  a graph  of  three  of  the  computed  modal  e fields 

z 

(eigenfunctions)  in  the  radiating  walls  of  the  pentagon  patch.  Since 

the  open-circuit  wall  impedance  condition  was  used,  the  mode  vectors  are 

entirely  real.  The  eigenvalues  and  eigenvectors  for  modes  27,  28  and  31 

2 2 2 

are  shown.  It  is  noted  that  eigenvalues  X = k00  = 0.1203  and  a„-,  = k„-. 

Zo  Zo  z/  11 

= 0.13605  are  closely  spaced,  indicating  the  possibility  of  a multiple- 
mode condition  being  produced  by  a single  feed. 

Fig.  4.21  illustrates  the  individual  magnetic  currents  in  each  wall 
for  modes  27  and  28.  The  net  effect  is  that  mode  28  is  vertically  polarized 
and  mode  27  is  horizontally  polarized.  The  relative  amplitudes  and  phases 
of  modes  27  and  28  would  be  determined  by  the  feedpoint  location  and  the 
frequency.  However,  these  preliminary  results  point  strongly  toward  a 
theoretical  prediction  of  circular  polarization  for  the  pentagon  patch,  if 
the  feedpoint  and  frequency  are  properly  chosen. 
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Fig.  4.17  Pentagon  microstrip  antenna  element  which 
produces  circular  polarization  (after 
Weinschel,  1975). 


4r  n702 


h of  the  modal  electric  fields  In. the  pentagon  radiating 
s as  a function  of  perimeter  position  (see  Fig.  4.19) 
modes  27,  28  and  31. 


PERIMETER  POSITION  (GRID  UNITS) 


0.S5x  4.?7 
= 3.7  3^ 


Fig.  4.21  Illustrating  individual  magnetic  currents  in  each 
wall  for  modes  27  and  28;  also  illustrating  vector 
sum  for  modes  27  and  28  - K and  K„  are  within 
1%  of  being  perpendicular. 


5.0  CONCLUSIONS 


This  report  has  presented  a method  of  analysis  which  can  be  applied 
to  microstrip  antennas.  The  primary  emphasis  to  date  has  been  on  the 
prediction  of  the  interior  fields,  since  this  is  a more  difficult  problem 
than  the  prediction  of  the  antenna  pattern. 

For  formulations  based  on  a cavity  model  with  open-circuit  walls 
(fields  satisfying  a Neumann  boundary  condition)  both  the  interior  mode 
vectors  and  eigenvalues  are  real.  However,  the  microstrip  patch  presents 
a wall  admittance  which  is  complex,  with  the  capacitive  susceptance 
usually  being  several  times  the  conductance.  Thus,  the  mode  vectors  are 
nearly  real  and  the  eigenvalues  are  nearly  real.  Solution  of  a trans- 
cendental equation  can  very  accurately  predict  the  resonant  frequencies 
once  the  appropriate  wall  admittance  is  known. 

The  problem  of  finding  the  wall  admittance  of  an  arbitrary  microstrip 
patch  shape  has  not  yet  been  solved  and  emphasis  should  be  placed  on  finding 
solutions.  However,  an  approximate  formula  has  been  found  which  uses  the 
conductance  of  a TEM  parallel-plate  waveguide  and  the  susceptance  of  an 
open-ended  microstrip  capacitor.  This  gives  good  agreement  between 
measured  and  calculated  input  impedances  for  the  rectangular  and  circular 
patches  considered.  The  input  resistance  at  the  resonant  frequency  is 
determined  by  the  wall  conductance,  whereas  the  resonant  frequency  is 
primarly  related  to  the  patch  dimensions  and  the  wall  susceptance. 

For  a patch  fed  by  a coaxial  probe,  the  input  impedance  to  the  patch 
proper  will  differ  from  the  measured  input  impedance  to  the  coaxial 
connector  by  a series  inductive  reactance  term,  given  by  eqn.  (4-24) . 

This  term  must  be  known  in  order  to  correct  measured  impedances  with  a 
probe  to  the  true  input  impedances  seen  by  a microstrip  transmission  line 
used  in  an  array  design. 

When  a corner-fed  nearly- square  patch  operates  in  a two-mode  config- 
uration with  two  closely-spaced  resonant  frequencies,  it  is  possible 
theoretically  to  obtain  circular  polarization  over  a 1%  bandwidth.  This 
multiple-mode  effect,  discussed  in  Sec.  4.1.3,  has  been  experiment- 
ally confirmed  at  PSL. 

The  method  of  finite  elements  has  been  applied  to  the  analysis  of  a 
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five-sided  patch  which,  when  properly  excited,  can  produce  circular 
polarization.  This  analysis,  which  is  in  a preliminary  stage,  confirms 
the  existence  of  two  spatially-orthogonal  exterior  fields  at  broadside 
and  at  two  closely-spaced  resonant  frequencies.  The  feedpoint  problem 
has  not  yet  been  solved  for  this  geometry. 

Both  the  classical  expansion  analysis  and  the  numerical  analysis 
approaches  will  be  continued.  The  finite  elements  program  (MICRO) 
will  be  expanded  to  include  (1)  sparse  matrix  algorithms  to  decrease 
computer  storage  requirements  and  (2)  subspace  iteration  to  decrease 
execution  time. 
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